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Abstract 

A  new  planning  procedure  is  proposed  for  economies  with  free  disposal 
in  production,  and  monotonic  social  welfare  functions.  No  convexity  of 
any  kind  is  required.  The  procedure  takes  the  novel  approach  of  searching 
for  an  optimum  in  the  region  of  the  unknown  production  set  that  the 
CPB  knows  the  least  about.  This  is  different  from  the  typical  approach  of 
searching  in  the  region  of  the  current  estimate  of  the  production  set  that  is 
most  preferred.  This  allows  an  upper  bound  to  be  calculated  on  the  speed 
of  convergence,  and  makes  it  possible  to  estimate  the  distance  between 
the  current  tentative  plan  and  the  optimal  plans.  This  is  an  essential 
feature  if  the  procedure  is  likely  to  be  stopped  before  it  converges,  and 
is  absent  in  existing  mechanisms.  The  procedure  is  also  able  to  identify 
and  discard  information  about  regions  of  the  production  set  that  cannot 
contain  optimal  plans. 


1.  Introduction 

With  the  changes  taking  place  in  the  Eastern  block,  one  might  be  led  to  wonder 
about  the  continued  relevance  of  economic  planning  mechanisms.  Most  governments 
are  turning,  at  least  in  part,  to  market  based  solutions  to  the  problem  of  economic 
development.  We  know  the  price  mechanism  is  well  suited  to  finding  efficient  al- 
locations in  convex  market  economies.  This  begs  the  legitimate  question  of  why  a 
planning  mechanism  should  be  used  at  all  in  such  economies.  Indeed,  some  plan- 
ning mechanisms,  Weitzman  (1970)  for  example,  are  highly  suggestive  of  a  market 
tatonement  process. 

Unfortunately,  convexity  is  an  unreasonably  strong  restriction  in  many  situa- 
tions. In  Public  Finance,  we  know  that  externalities  lead  to  fundamental  noncon- 
vexities  in  the  production  sets  (see  Starrett  (1972)).  In  International  Trade,  there 
are  concerns  about  increasing  returns  and  learning  by  doing.  In  Industrial  Orga- 
nization, nonconvexities  arise  from  natural  monopoly  and  other  forms  of  market 
power.  Indivisibilities  in  large  public  projects  can  also  lead  to  nonconvex  produc- 
tion sets.  In  such  cases,  we  have  no  assurance  that  the  competitive  equilibrium  will 
even  exist.  Given  this  failure  of  the  market,  an  alternative  means  must  be  employed 
to  find  socially  optimal  allocations. 

This  paper  describes  an  iterative  planning  procedure  which  may  be  used  in 
a  broad  class  of  economies.  We  require  only  that  each  firm  in  the  economy  has 
a  closed  production  set  and  produces  under  conditions  of  free  disposal,  and  the 
Central  Planning  Board  (CPB)  has  a  continuous  and  monotonic  Social  Welfare 
Function  (SWF).  We  make  no  convexity  assumptions  of  any  kind. 

Following  the  tradition  of  Malinvaud  (1967),  Weitzman  (1970),  and  especially 
Cremer  (1977),  we  frame  the  problem  as  a  CPB  trying  to  maximize  a  SWF  over  an 
unknown  feasible  set.  To  do  this,  the  CPB  asks  a  series  of  simple  questions  of  the 
firms  in  the  economy  in  order  to  generate  successively  better  approximations  of  the 


feasible  set. 

For  the  one  firm,  two  good  case,  the  procedure  is  roughly  as  follows  (It  may 
help  to  refer  to  figure  1.)  The  CPB  starts  off  by  choosing  a  point  that  overesti- 
mates and  another  point  that  underestimates  a  part  of  the  Production  Possibilities 
Frontier  (PPF).  Ideally  this  part  of  the  PPF  contains  an  optimal  plan,  but  this  is 
not  required.  We  can  imagine  that  the  overestimate  and  the  underestimate  are  the 
corners  of  a  box  around  part  of  the  PPF.  Call  the  overestimate  the  "best  point"  of 
the  box  and  the  underestimate  the  "worst  point"  of  the  box.1  We  will  show  that 
the  procedure  converges  to  the  most  preferred  point  on  the  PPF  in  this  initial  box.2 

Boxes  have  two  edges  radiating  out  of  the  best  point  (and  of  course  two  sym- 
metric edges  coming  from  the  worst  point.)  This  initial  box  is  iteratively  subdivided 
into  a  set  of  smaller  boxes  that  surround  the  PPF  in  following  way.  Each  iteration, 
the  CPB  finds  the  box  whose  longest  edge  is  longer  than  every  other  box's  longest 
edge  (at  the  first  iteration,  there  is  only  one  box  and  so  the  choice  is  trivial.)  We 
take  this  "biggest"  box  to  be  the  region  the  CPB  knows  the  least  about.  The  CPB 
bisects  the  longest  edge  of  the  biggest  box  with  a  perpendicular  plane  and  asks  the 
firm  to  report  a  point  in  the  intersection  of  the  PPF  and  this  plane.  If  such  a  point 
exists,  the  box  is  divided  in  half  along  its  longest  edge.  Note  that  this  also  adds  a 
new  best  point  and  worst  point.  If  no  such  feasible  point  exists,  it  must  be  that  the 
PPF  either  lies  entirely  above  or  entirely  below  the  perpendicular  plane  within  the 
box.  In  this  case,  the  biggest  box  is  divided  in  half  and  whichever  of  the  two  boxes 
does  not  contain  any  part  of  the  PPF  is  discarded.  Thus,  each  iteration,  the  box 
that  is  currently  the  largest  is  divided  in  half  along  its  largest  edge. 


We  call  the  overestimate  the  best  point  since  by  monotonicity  it  is  the  most  preferred  point  in  the 
box.  The  reverse  is  true  for  the  worst  point. 

Thus,  by  choosing  a  larger  initial  box,  the  CPB  makes  it  likely  that  the  procedure  will  converge 
to  a  better  plan.  However,  this  also  slows  down  convergence,  as  a  larger  part  of  the  PPF  has  to 
be  considered. 


The  procedure  uses  the  information  acquired  above  as  follows.  The  points 
on  the  PPF  that  the  firm  reports  allow  the  CPB  to  find  an  underestimate  (called 
"lower  estimate"  in  this  paper)  of  the  production  set.  By  free  disposal,  the  CPB 
knows  that  everything  that  lies  below  a  point  that  is  feasible  is  also  feasible.  Thus, 
the  comprehensive  hull  3  of  all  of  the  known  feasible  points  may  be  taken  as  the 
lower  estimate.  The  set  of  best  points  described  previously  is  constructed  such  that 
taking  their  comprehensive  hull  gives  a  set  that  contains  the  PPF  and  thus  is  an 
overestimate  (called  "upper  estimate"  in  this  paper)  of  the  production  set.  The 
most  preferred  point  on  the  lower  estimate  at  any  given  iteration  is  a  feasible  point 
which  is  taken  as  the  current  tentative  plan.  Note  that  since  this  lower  estimate 
expands  each  iteration,  the  set  of  tentative  plans  is  monotonically  increasing  in  the 
preference  ordering.  The  procedure  has  one  last  feature  which  reduces  its  informa- 
tional requirements.  If  there  is  a  best  point  that  is  inferior  under  that  SWF  to  a 
tentative  plan,  than  that  best  point  clearly  cannot  lie  above  an  optimal  plan.  Thus 
the  CPB  may  discard  boxes  dominated  by  such  inferior  best  points.4 

In  summery,  the  algorithm  proceeds  as  follows.  Each  iteration,  the  CPB  finds 
the  most  preferred  point  on  the  current  lower  estimate  and  names  this  the  tentative 
plan.  It  then  discards  all  boxes  whose  best  point  is  inferior  to  the  tentative  plan. 
Next,  the  CPB  finds  the  "biggest"  box  remaining  and  divides  it  in  half  along  its 
longest  edge.  If  both  of  the  two  new  boxes  contain  parts  of  PPF,  then  both  are 
retained,  and  the  firm  reports  any  point  in  the  intersection  of  the  PPF  and  the 
plane  that  bisects  and  is  perpendicular  to  the  longest  edge.  Appropriate  new  points 
are  added  to  the  set  of  best,  worst  and  feasible  points.  Otherwise,  the  box  that  does 
not  contain  a  part  of  the  PPF  is  discarded  and  the  set  of  best  and  worst  points  is 
updated  appropriately. 


See  the  next  section  for  a  precise  definition  of  comprehensive  hull. 

This  identification  of  irrelevant  information  is  a  novel  feature  of  this  mechanism,  as  far  as  we  know. 


Most  planning  procedures  ask  firms  to  refine  the  CPB's  information  in  the 
region  of  the  current  estimate  of  the  production  set  that  is  most  preferred  under  the 
SWF.  However,  in  a  nonconvex  environment,  there  is  no  particular  reason  to  believe 
that  the  optimal  plan  is  close  to  the  currently  most  preferred  region.  Indeed,  the 
indifference  curve  through  the  optimal  plans  may  come  very  close  to  many,  widely 
dispersed  areas  of  the  PPF  if  neither  is  convex.  In  the  mechanism  described  in 
this  paper,  the  CPB  asks  each  firm  to  refine  the  board's  information  in  the  region 
it  knows  the  least  about.  This  is  a  novel  approach,  as  far  as  we  know,  and  yields 
several  advantages. 

In  general,  most  mechanisms  cannot  be  expected  to  converge  to  an  optimal 
plan  in  a  finite  number  of  iterations.  Even  if  this  were  not  the  case,  Bennett  (1985) 
points  out  that  time  shortages  may  require  that  the  procedure  be  stopped  before 
convergence,  and  the  current  tentative  plan  taken  as  the  final  approximation  of 
an  optimal  plan.  Given  this  real  world  constraint,  it  is  essential  that  CPB  know 
how  quickly  the  various  planning  algorithms  start  suggesting  tentative  plans  that 
are  satisfactorily  close  to  an  optimal  plan.  Without  this  information,  there  is  no 
reasonable  way  to  choose  which  algorithm  to  use. 

This  problem  has  led  several  authors  to  investigate  the  question  of  speed  of 
convergence.  One  approach  is  to  show  that  one  procedure  is  always  faster  than  an- 
other. For  example,  Clark  (1989)  compares  the  procedures  of  Malinvaud  (1967)  and 
Chander  (1978).  He  shows  that  Malinvaud's  is  at  least  as  fast  as  Chander's  in  a  cer- 
tain class  of  economies,  although  Malinvaud's  procedure  has  a  higher  informational 
requirement.  Unfortunately,  it  is  very  difficult  to  generalize  this  type  of  comparison 
to  other  procedures.  It  would  be  quite  surprising  if  it  were  possible  to  rank  all  the 
known  procedures  on  the  basis  of  speed,  given  the  diversity  of  approaches  taken  in 
this  literature.  In  any  event,  such  theorems  still  leave  the  CPB  ignorant  about  the 
actual  speed  of  convergence  of  any  given  mechanism. 

An  alternative  approach  is  to  conduct  computer  simulations  in  order  to  test 


which  procedure  is  faster.  Clark  gives  some  very  interesting  results  of  such  an  ex- 
periment with  the  Malinvaud  and  Chander  procedures.  However,  as  Chander  and 
Kundu  (1991)  point  out,  the  results  are  highly  dependent  on  the  class  of  examples 
considered.  To  be  able  to  make  useful  assertions  like  "procedure  a  converges  faster 
than  procedure  /?,  x  percent  of  the  time,"  we  would  need  to  have  a  measure  over 
the  class  of  reasonable  economies.  Denning  this  class,  much  less  a  sensible  mea- 
sure, seems  an  impossible  task,  especially  if  we  do  not  restrict  attention  to  convex 
economies.  So  while  these  results  are  suggestive,  and  interesting,  it  would  be  overly 
hasty  to  choose  a  procedure  based  on  simulations  alone. 

We  take  a  different,  and  we  argue,  much  more  complete  approach  in  this  paper. 
It  is  easy  to  find  an  upper  bound  on  the  longest  edge  of  the  biggest  box  at  any  given 
iteration.  Since  this  tells  us  the  maximal  size  of  the  boxes  that  surround  the  PPF, 
we  can  use  this  to  find  an  upper  bound  on  the  speed  of  convergence  in  several 
different  ways.  This  bound  may  then  be  employed  to  make  direct  comparisons  with 
other  procedures  for  which  a  bound  on  speed  of  convergence  is  known.  In  addition, 
the  actual  degree  of  ignorance  (size  of  the  largest  box)  may  be  calculated  exactly  as 
the  algorithm  progresses.  This  makes  it  possible  to  set  up  informed  stopping  rules 
of  the  form:  stop  searching  for  an  optimum  if  the  current  tentative  plan  is  within  a 
given  distance  of  an  optimal  plan  under  some  measure  (eg.  within  a  certain  utility 
value,  or  within  a  certain  e-ball  in  the  goods  space). 

Thus,  the  major  advantage  of  this  mechanism  is  that  the  novel  method  of  infor- 
mation acquisition  makes  it  possible  to  directly  estimate  the  speed  of  convergence, 
and  to  set  up  sensible  stopping  rules.  It  is  not  clear  how  one  would  do  this  for 
existing  procedures  denned  for  nonconvex  economies  such  as  Henry  and  Zylberberg 
(1978)  or  Cremer  (1977,  and  1983). 


2.  Definitions  and  assumptions 

This  section  lays  out  the  basic  definitions  and  assumptions  that  are  used  in 
what  follows.  We  consider  a  CPB  with  a  complete  and  transitive  social  preference 
ordering,  ^,  over  a  closed  consumption  set  with  m  goods:  X  C  Rm .  The  CPB  is 
endowed  with  u>  G  3?m  at  the  outset.5  We  assume: 

Al)  The  preference  relation  y  is  continuous. 
A2)  If  x  >  x'  then  x  >-  x' .6 

(Preferences  are  strongly  monotonic.) 

The  weak  upper  and  lower  contour  sets  of  this  relation  for  any  subset  Z  C  Rm 
are  denoted  by  U(Z)  and  L(Z)  respectively.  Formally: 

U(Z)  =  {x  G  X  |  x  y  x  for  some  x'  G  Z} 

and, 

L(Z)  =  {x  G  X  |  x  ■<  x'  for  some  x   G  Z)  . 

The  boundary  of  a  set  is  denoted  by  "5".  Under  these  two  assumptions,  dL(Z)  is 
the  indifference  surface  through  the  infimum  of  Z  under  the  preference  ordering. 
These  assumptions  also  imply  that  if  x  G  X,  and  x'  >  x  then  x  €  X. 

There  are  k  firms  in  the  economy  with  production  sets  denoted  by  Ye  for 
£  =  1, . . . ,  k.  We  make  two  assumptions  on  firms'  production  sets: 

Bl)  If  y  G  Yl ,  then  y'  <  y  implies  y'  G  Yl  for  £  =  1, . . . ,  k. 

(Free  disposal  in  production,  or  identically,  comprehensiveness  of  Yt .) 


5  It  is  possible  to  reformulate  the  problem  such  that  the  firms  own  the  initial  endowments  without 
changing  the  procedure  in  any  substantial  way.  However,  this  is  a  planning  problem,  so  there  does 
not  seem  to  be  any  gain  from  doing  so. 

The  three  types  of  vector  inequalities  are  3>,  >,  and  >  . 


B2)   Yl  is  closed. 


We  dispense  with  all  convexity  assumptions,  and  depend  on  monotonicity  and 

free  disposal  instead.  Without  convexity,  closedness  of  each  firms'  feasible  set  is  not 

enough  to  conclude  that  the  set  of  feasible  allocations  for  the  economy  is  also  closed 

in  general.   In  this  algorithm,  however,  the  search  for  an  optimal  production  plan 

is  limited  to  a  compact  subset  of  the  each  firms'  production  set.  Since  the  sum  of 

these  compact  subsets  of  Yl  is  the  only  part  of  the  global  production  set  that  is 

considered,  the  relevant  part  of  the  global  production  set  is  compact. 

The  usual  notion  of  set  summation  is  used: 
t 
^Zl  =  {x  e$lm  \z  =  zl  +---  +  Z1  <mdzl  £  Zx  fort  =  l,...,t}. 
t=i 

Each  global  object,  z,  (a  global  feasible  production  plan,  for  example)  is  the  sum 
of  one  element  from  each  of  the  Ze's.  That  is,  z  =  z1  +  . . .  +  zk  where  zl  €  Zl  for 
I  =  1, . . . ,  k.  We  denote  by  z  this  set  of  k  elements  that  sum  to  z.  For  example, 
if  we  use  this  notation  to  indicate  the  decomposition  of  a  set  of  vectors,  Z,  this 
becomes: 


r  *i  ^ 


Z  =  <   Zi 


»  za  t 


/      1  2 


>  =  <  zj 


,*-!         -rk  \ 


Z*"1        Zk 
I  i 


k  z1     z2 

<      s  s 


Zs      X        Zs    , 


\  e% 


mx  kxs 


where  zj  £  Ze  for  all  j  and  all  I  =  1, . . . ,  fc,  and  £^=1  ze-  =  Zj. 

Since  we  will  not  consider  the  case  of  production  externalities,7   the  set  of 
feasible  plans  for  the  economy  is: 


Y={yEVT 


1=1  ) 


This  generalization  is  easily  accomplished  without  any  change  in  the  model  except  the  definition 
of  the  global  feasible  set.  However,  the  notational  cost  (firm  i's  observation  of  firm  j's  output  of 
good  k)  would  be  heavy  in  an  already  notationally  dense  paper. 


This  simply  says  that  a  plan  is  feasible  for  the  economy  if  it  is  the  sum  of  the 
endowment  and  a  global  production  plan  which  can  be  decentralized  in  such  a 
way  that  each  firm  is  asked  to  produce  a  feasible  net-output  vector.  We  will  use  the 
convention  that  superscripts  refer  to  the  firms,  and  omission  of  superscripts  indicate 
a  global  object.  Thus  ye  6  Yl  is  a  production  plan  for  firm  £,  while  y  E  Y  is  a 
production  plan  for  the  economy. 

The  procedure  depends  heavily  on  the  notions  of  taking  comprehensive  hulls 
of  various  sets.  For  any  set  Z  C  9£m,  define  the  comprehensive  hull  of  Z  as  follows: 

ch(Z)  e  {x  G  r  |  3  2  €  Z  s.t.  x  <  z) 

This  is  the  set  of  vectors  in  9ftm  that  are  weakly  dominated  by  some  vector  in  Z.  We 
will  also  have  need  of  the  notion  of  an  inverse  comprehensive  hull  of  Z  in  explaining 
the  mechanism: 

ich(Z)  =  {xe^n\^zeZ  s.t.  x  >  z} 

This  is  the  set  of  vectors  in  9£m  that  do  not  strictly  dominate  some  vector  in  Z. 
Note  that  both  the  comprehensive  hull  and  inverse  comprehensive  hull  of  a  finite 
collection  of  vectors  is  closed. 

Finally,  consider  the  following  notion  of  minimal  "distance":8 

u{A,B)=  infinf   \\  x  -  y  \\ 

y(zB  iG-A 
where  A  and  B  are  sets  in  9Rm,  and  ||  •  ||  denotes  the  Euclidean  norm,  v  gives  the 
infimum  of  length  of  the  smallest  gap  between  two  sets.  Thus,  if  A  and  B  intersect, 
then  v(A,  B)  =  0.  A  useful  property  of  v  is  given  in  Lemma  1. 

Lemma  1.  If  A'  C  A  C  &m  and  B'  C  B  C  ftm,  are  all  bounded  then,  v(A',B')  > 
HA,B). 

Proof/ 


Note  that  u  is  not  a  measure  of  distance  in  a  mathematical  sense  since  it  does  not  satisfy   the 
triangle  inequality. 


By  the  definition  of  i/,  and  of  compact  sets,  there  must  be  two  points,  a'  E 
closure(A'),  and  V  €  closure(B')  such  that  v(A',B')  >  v(a',b').  But  then,  a1  £ 
closure^)  and  V  €  closure(B).  So  u(A',B')  >  u(a',V)  >  v(A,B). 


3.  The  Algorithm 

The  algorithm  is  based  on  the  observation  that  it  is  possible  to  construct  a 
lower  bound  as  well  as  an  upper  bound  on  the  production  set  of  each  firm  using  the 
information  acquired  through  the  Cremer  procedure.  These  estimates  can  then  be 
used  jointly  to  narrow  the  area  of  search  with  each  iteration.  Basically,  this  is  done 
by  removing  from  consideration  all  parts  of  the  overestimate  of  the  production  set 
that  are  inferior  to  some  part  of  the  underestimate.  In  the  Cremer  procedure,  and  in 
most  other  iterative  planning  procedures,  the  entire  feasible  set  must  be  considered 
as  potentially  containing  an  optimum  unless  the  procedure  actually  converges. 

Cremer  constructs  an  overestimate  of  the  production  set  of  each  firm  by  using 
free  disposal  to  conclude  that  points  that  strictly  dominate  elements  of  the  PPF 
cannot  themselves  be  feasible.  Formally,  he  constructs  an  Upper  Bound  on  the 
Production  Set  of  the  ith  firm  at  iteration  n  thus: 

VBPFen  =  zch(Xen)f]ch(beQ) 

where  X*  is  the  set  of  points  known  to  be  on  the  d.th  firm's  being  PPF  at  the  nth 
iteration,  and  6q  is  a  vector  that  dominates  at  least  one  optimal  production  plan. 
Note  that  UBPF„  is  a  closed  and  comprehensive  set  as  the  intersection  of  two  such 
sets. 


A 


Xn  —  {Xj  ,  x2,  x3,  x4,  x5j 

Figure  1 

Free  disposal  can  obviously  be  applied  in  the  other  direction  as  well.  If  a  point, 
x,  is  known  to  be  feasible,  then  all  points  that  x  dominates  must  also  be  feasible. 
The  CPB  may  therefore  take  the  comprehensive  hull  of  all  known  feasible  points  as 
a  lower  bound  on  the  production  set.  Thus,  we  define  the  Lower  Estimate  of  the 
Hh  firm's  production  set  at  iteration  n  as: 

LEen  =  ch(Xen). 


Figure  1  illustrates  these  upper  and  lower  bounds  on  a  firm's  production  set. 

10 


Notice  that  the  set  theoretical  difference  between  them  is  a  union  of  "cubes".9  This 
will  also  be  true  in  higher  dimension.  Each  cube  is  uniquely  characterized  by  a 
point  6,  which  dominates  all  other  points  in  the  cube  (called  the  best  point),  and  a 
point  w,  which  is  dominated  by  all  other  points  in  the  cube  (called  the  worst  point) 
via  the  correspondence: 

C(b,w)  =  {xe^m\b>x>w}. 

Observe  that  Cremer's  upper  bound  on  the  production  set  can  also  be  found  by 
taking  the  comprehensive  hull  of  all  the  best  points.  Since  the  set  of  best  points 
is  explicitly  constructed  in  the  algorithm  described  in  this  paper,  it  is  natural  that 
the  upper  bound  on  the  production  set,  called  the  Upper  Estimate  here,  should  be 
defined  accordingly: 

UEn  =  ch(Ben). 

The  algorithm  itself  involves  the  iterative  updating  of  four  sets  of  quantity 
vectors  and  the  naming  of  a  tentative  production  plan,  pen,  for  each  firm,  for 
each  iteration.  The  first  of  these  is  C£,  a  set  of  ordered  pairs  {(6(,u>(),  ..., 
(&*, tuf), . . . ,  }  =  {c\ ,,...,  cf, ...,  }  such  that  when  these  pairs  are  used  to  form 
cubes,  their  union  equals  the  set  theoretical  difference  between  U Eln  and  LE„.  The 
next  two  are  Bln  and  W%,  the  sets  of  best  and  worst  points  in  Cln  (that  is,  the 
set  of  all  first  and  second  elements,  respectively,  in  the  set  of  ordered  pairs  C„). 
It  is  notationally  convenient  to  construct  these  sets  separately  even  though  all  the 
information  contained  in  them  is  also  contained  in  Cln.  Last  is  X^,  which  is  a  set 
of  points  known  to  lie  exactly  on  firms'  PPF's. 

All  of  the  sets  above  are  constructed  individually  for  each  firm.  These  are 
aggregated  to  form  overestimates  and  underestimates  of  the  global  PPF  as  follows: 


The  use  of  the  word  "cube"  is  not  meant  to  imply  that  these  objects  have  equal  sides,  or  that 
they  are  three  dimensional.  We  use  "cube"  as  a  substitute  for  the  more  accurate,  but  awkward 
"hyper-  rectangle" . 
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*-  =  £« 

e=i 

UEn  =  ch(Bn) 

LEn  =  ch(Xn). 

Thus,  the  global  upper  and  lower  estimates  are  of  the  set  of  feasible  global  net 
output  vectors,  not  of  the  production  set  itself.  The  set  of  global  feasible  allocations 
is  found  by  adding  the  u>  to  the  upper  and  lower  estimates.  We  follow  this  convention 
for  notational  convenience. 

The  CPB  must  make  an  initial  over  and  underestimate  of  the  production  set 
for  each  of  the  firms.  These  initial  estimates  may  be  as  broad  or  narrow  as  desired. 
We  show  that  the  algorithm  converges  to  the  best  plan  possible  within  these  limits. 
In  particular,  if  we  assume  that  the  CPB  knows  enough  to  be  able  to  over  and 
underestimate  a  part  of  the  production  set  that  contains  an  optimal  plan,  then  the 
algorithm  converges  to  the  true  optimum.  Cremer  makes  an  assumption  in  this 
spirit.  Our  approach  is  make  explicit  the  trade  off  between  speed  of  convergence 
and  the  possibility  of  finding  a  better  plan  that  choosing  the  boundaries  of  search 
entails. 

Formally,  for  each  firm,  the  CPB  chooses  initial  over  and  underestimates  of  the 
production  sets,  (&q,u>q),  such  that: 

V*=l,...,Jb,  beQ  (£Ye,we0  EY',  and&J  >>  we0. 

The  algorithm  is  then  initialized  as  follows: 

Cfs{cJ}s{(fi{,wX)}for/  =  lf...1*1 
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4  =  {we0}  for  I  =  1, . . .  ,  k  X0  =  {w0}  =  I  Y,  wo  | 

Po  =  Wo)  for  *  =  1, . .  • ,  *  Po  =  {w0}  =  <  ^  wlQ  \  , 

where  p$  1S  ^e  tentative  production  plan  for  firm  I  at  iteration  0,  and  po  is  the 
corresponding  global  tentative  plan. 

Neither  the  set  of  global  cubes  nor  the  set  of  global  worst  points  is  collected. 
This  is  because  the  global  best  points  are  constructed  in  such  a  way  that  their 
comprehensive  hull  gives  an  overestimate  of  the  production  set.  The  comprehensive 
hull  of  the  known  points  Xn,  on  the  other  hand,  gives  an  underestimate.  These  two 
estimates  are  all  that  we  need  at  the  global  level.  Cubes  are  useful  only  in  that 
their  size  may  be  used  as  a  measure  of  the  CPB's  ignorance  about  a  particular  firm's 
production  possibilities.  The  CPB  will  end  up  asking  the  firms  questions  about  the 
part  of  the  PPF  contained  in  the  "biggest"  cube.  Since  no  analogous  question  is 
ever  asked  at  the  global  level,  global  cubes,  and  consequently,  global  worst  points 
are  not  needed. 

We  will  show  that  the  algorithm  converges  to  the  best  plan  within  the  search 
boundaries  defined  by  above.  Formally,  let  Y  :  3?  m  — >  3ftm  be  the  correspondence 
that  gives  the  subset  of  feasible  global  net  output  vectors  examined  when  the  CPB 
chooses  Co,  as  his  search  area: 

Y(c0)  =  Y(cl...,ck0)  =  L  =  J2y^    V£=l,...,A:,y/e{c(6{,u;J)f|r/}|. 

Then  let  X*  :  3?2m  x  — >  $Rm  be  the  correspondence  that  gives  the  set  of  optimal, 
global  net  output  vectors  in  Y(co): 

X*(c0)  =  {ye  Y(c0)  \y  +  uyy'  +  uVy'e  Y(c0)} . 
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Again,  these  correspondences  give  feasible  and  optimal  net  output  vectors,  not 
allocations.  We  are  finally  ready  to  define  the  algorithm.  Each  iteration  is  broken 
down  into  five  steps.  Figure  2  gives  an  illustration  for  the  one  firm  case. 

Step  1)  The  first  step  for  any  given  iteration,  n,  is  to  take  the  most  preferred 
point  on  the  global  lower  estimate  for  the  previous  iteration,  LEn-\,  as  the  global 
tentative  plan  pn.  Since  LEn-\  =  ch(Xn-i),  and  the  CPB's  preferences  are  mono- 
tonic,  this  maximization  must  take  place  at  some  element  of  Xn-\.  Each  element 
x  £  Xn-\  has  a  known  decomposition  x  =  (x1, . . .  ,x*).  Thus,  pn  may  be  decom- 
posed into  pn,  which  in  turn  is  taken  as  a  specification  of  a  tentative  production 
plan  for  each  firm.  Formally: 

1)  Some  pn  £  {x  £  LEn-\  |  i+w  >:  t/+w  V  y  G  LEn-\)  is  chosen  and  pn  is 
declared  to  be  the  tentative  plan  for  the  economy. 

Step  2)  The  second  step  is  to  find  all  elements  of  -0n-i  that  are  strictly  inferior  to 
pn,  and  discard  them  to  form  the  new  set  Bn.  We  do  this  because  any  element  of 
the  overestimate  that  is  inferior  to  a  feasible  point  pn  certainly  cannot  lie  above  an 
optimal  plan.  Such  points  may  therefore  be  safely  removed  from  future  considera- 
tion. Likewise,  elements  that  are  known  to  be  feasible,  but  which  are  nevertheless 
dominated  by  a  discarded  element  of  J5n,  may  be  thrown  away.  Formally: 

2)  Bn  =  {b  £  Bn-X  |  pn  +  u  ■<  b  +  u} 

Xn  =  <  x  6  Xn-\  |  b  >  x  for  some  b  E  Bn  > 

Step  3)  Having  found  and  discarded  the  irrelevant  elements  of  the  global  sets  Bn-\ , 
and  Xn-\,  the  next  step  is  to  find  the  irrelevant  elements  of  the  sets  collected  from 
each  of  the  firms.  Suppose  that  a  best  point  for  a  firm,  be  £  Bin_l,  is  only  used 
to  create  global  best  points  which  are  known  to  be  inferior  to  a  feasible  point. 
Then  clearly,  the  firm  will  never  be  called  upon  to  produce  be ,  or  any  point  that 
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dU(X*(co)  +  w) 


dU(X*(co)+tj) 


d£(Pn  +  U 


dL(Pn  + 


c~ 


2a)  Step  1:  xx  is  the  most  prefered 
point  on  LEn  and  is  declared  the 
tenative  plan,  pn. 


dU(X*(co}+«>) 


2b)  Steps  2  and  3:  b3  is  found  to 
be  inferior  to  pn  and  so  63,  u;3, 
and  c3  are  discarded. 


dU(X*{co)  +  u) 


w2 

The  point 
^0     reported  by 
!  X         the  firm 


Y<MPn  +  Cj) 


\^dL(vn  +  c 


D 


2c)  Step  4:  The  xth  edge  of  the  second  2d)  Step  5:  c2  is  divided  in  half  aW 

cube  is  found  to  be  the  longest.  the  x-axis  and  replaced  with  the  half 

Case  1  holds,  and  the  firm  reports  new  cubes,  and  appropriate  best  and 

a  point  in  the  in  the  intersection  of  worst  points  added  to  Bn  and  W 
n{c2)  and  the  PPF 


Figure  2:  The  new  procedure 
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b£  dominates.  This  is  identical  to  saying  that  if  be  is  only  used  to  create  global 
points  that  are  not  in  the  set  Bn,  then  it  may  be  discarded.  Furthermore,  all  points 
of  the  sets  W^_x  and  C^l_1  that  are  associated  with  discarded  best  points,  and  all 
elements  of  X^l_1  that  are  dominated  by  discarded  points  of  Bfl_1  may  be  discarded. 
Formally: 

3)  (*!,...,**)  s(£n),10 

and  for  each  £  =  1, . . . ,  k: 

W%  =  {we  £  Wln_x  |  3  (6'X)  =  cl  6  Cln_x  where  be  (E  Ben} 

X£n  =  {xe  e  Xln_x  \be>xe  for  some  be  6  Ben}. 

Step  4)  Having  discarded  irrelevant  information  about  each  of  the  firms'  production 
sets,  the  next  step  is  to  gain  more  information  about  the  parts  that  remain.  In  this 
algorithm,  the  CPB  asks  for  information  about  the  region  he  knows  the  least  about. 
Notice  that  each  cube  has  m  edges  radiating  out  from  its  best  point.  The  CPB  is 
interested  in  the  cube  whose  longest  edge  is  longer  than  every  other  cube's  longest 
edge.  Let  the  function  ILE  :  3£m  x  5Rm  — >  {1, . . . ,  m}  give  the  Index  of  the  Longest 
Edge  of  a  cube  where  ties  are  broken  by  taking  the  lowest  index: 

ILE{ce)=  {j  G(l,...,m)  |  Vi  =  l,...,m  &>'•'  -  wJ'e  >  b4'1  -  wl'e, 

and  j  <  i  Vij'  =  l,...,m  s.t.  b3^  -  w3'e  =  bl'e  -  wl,i}  . 

In  all  cases  the  first  superscript  refers  to  a  component  of  a  vector  while  the  second 
refers  to  the  firm.  Let  E„*  denote  the  length  of  the  longest  edge  any  cube  in  C„: 

Een*  =  lb3  -  wJ  |  b3  -w3  >bl  -wlVc,ce  Cen  and  V i,j  =  1, . . .  ,m}  . 


Any  given  element  6  may  be  summed  with  many  different  permutations  of  best  points  from  other 
firms  to  form  the  set  of  global  best  points.  Thus,  when  the  set  of  global  best  points  is  decomposed, 
there  may  be  a  great  deal  of  repetition  in  the  list  of  best  points  for  each  firm. 
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Note  that  if  E„*  is  the  length  of  the  longest  edge  of  any  cube  in  C^,  it  is  also  an 

upper  bound  on  the  length  of  the  longest  edge  of  any  cube  in  C„,  since  no  cube 

grows  in  size  as  C„  is  converted  into  C„.   Correspondingly,  E„  is  an  upper  bound 

on  the  length  of  the  longest  edge  of  any  global  cube: 

Jt 
E*  =  V  Ee*. 

£=1 

Let  C£*  be  the  set  of  cubes  in  C„  whose  longest  edge  is  at  least  as  long  Ee*: 

Cen*  =  lc£Cln\V  -w>  =  El*  and  j  =  ILE(c)}  . 

For  any  cube  c,  consider  the  hyperplane  which  bisects  the  cube,  and  is  perpen- 
dicular to  the  longest  edge.  We  denote  the  correspondence  that  gives  this  hyperplane 
by  h  :  ft2m  —  3tm: 

h{ce)  =  ix  €  3?m  |  fori  =  ILE(ce),xJ  =  ^'   +  wh  ^ 

In  the  algorithm,  for  each  firm,  an  arbitrary  element  ce  6  C^*,  is  chosen  by  the 
CPB,  and  the  firm  is  asked  to  report  some  point  on  the  PPF  within  the  cube  ce ', 
and  on  h(ce).  That  is  to  say,  the  firm  is  asked  to  report  some  point  in  the  set: 

PPF(ce)  =  {h(ce)ndYenC{be,we)}  . 

If  such  a  point  exists,  it  is  added  to  the  set  X„  to  create  X„-  The  set,  PPF(ce), 
however,  may  be  empty.  In  this  case,  the  PPF  must  be  either  completely  above  or 
completely  below  the  hyperplane  within  the  cube.  This  is  because  of  the  assumption 
that  Yl  is  closed,  which  implies  that  the  PPF  is  continuous.  If  either  one  of  the 
above  is  true,  the  firm  reports  this  instead  of  a  point.  In  both  cases,  X^  is  set  equal 
to  X„.  It  is  now  possible  to  conclude  that  the  procedure  is  well  defined  since  this  is 
the  only  question  that  the  firms  are  ever  asked,  and  cases  (a),  (b),  and  (c),  below, 
are  exhaustive.  Formally: 
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4)  For  each  I  =  1, . . . ,  &,  ce  £  C„*  is  chosen,  and  the  firm  is  asked  to  report  any 
xe  £  PPF(ce).  The  firm  makes  one  of  three  responses: 

(a)  PPF(ce)  7^  0,  and  the  firm  reports  some  element  of  the  set,  x£ . 

(b)  PPF(c£)  =  ®<mdx<yVx€  C{be,we)  n  h(ce)  and  y  6  dYe  n  C(6',  w')- 

(c)  PPF(ce)  =  0andi>yVxe  C(6',  u/)  n  /i(c€)  and  y  £  dYe  H  C(be,we). 
In  case  (a),  Xen  =  lxenUxe\. 

In   case  (b)  and  (c),  X„  =  X„. 

Step  5)  The  fifth  and  final  step  of  the  algorithm,  which  is  illustrated  in  figure  3,  is 
to  update  the  sets  B^l_l,  W^_l,  and  C„_j.  This  was  partially  accomplished  in  step 
3  when  these  sets  were  turned  into  B„,  etc.  The  updating  is  completed  in  different 
ways  depending  on  what  the  firm  reports  at  step  4.  The  first  possibility  is  that  the 
set  PPF(ce)  is  not  empty,  and  the  firm  reports  some  element  xe  to  the  CPB.  In 
this  case,  the  cube  c£  is  divided  in  half  along  the  j  axis,  and  ce  is  replaced  with 
the  two  new  cubes.  Appropriate  best  and  worst  points  are  added  to  create  B„  and 
W*.  The  second  possibility  is  that  PPF(ce)  is  empty  and  h(c£)  C\C(be ,  we)  is  below 
dYe  (lC(be ,w*).  In  this  case  c  is  truncated  by  moving  its  worst  point  halfway  up 
the  longest  edge  of  the  cube.  The  last  possibility  is  that  PPF(ce)  is  empty  and 
h(ce)C\C(be ,  w£)  is  above  dYeC\C(be ,  we).  Here,  the  opposite  is  done.  The  best  point 
of  the  cube  is  moved  halfway  down  the  longest  edge.  Formally: 

5)  For  each  £=  1,  ...,&,  if  case  (a)  obtains: 
Bi  =  {(Bi\  b<)Ub<aUb<,} 

where  beQ  =  be , 

and  b^  =  (b1'1, . . . ,  &J-1'',  v-ly'\  V*1'*, . . . ,  6m'<); 

where  wla  =  (wl>e,  w2>£, . . .  w>-M   *,,<+*i,l>tpi+Mi . . .  ,  wm^ 
and  wi  =  we; 

ci^{(ci\  c')u4u4}. 
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D 


7YTbe  =  bi  y 


Wp  =  W 


PPF(ce) 
xen,  the  point  the  firm  reports. 

3a)  In  this  example,  assume  that  E„*  =  (bl,e  —  wl,e).  Then  the  cube  ce  is  to  be 
divided  according  to  case  (a)  since  PPF(c()  ^  0.  As  can  be  seen  in  the  figure,  c 
is  divided  exactly  in  half  along  the  first  edge  creating  two  new  cubes  ceQ  and  cl, 
whose  union  equals  c  . 


h(c<) 


0 


3b)  Here.  PPF(ce)  =  0  and  dYe  D  ce  is    3c)  Here,  PPF{ce)  =  0  and  dYe  n  c(  is 
below  h(ce).  Case  (b)  obtains,  and  we      above  h(ce).  Case  (c)  obtains,  and  be 
is  moved  up  to  w^  and  the  dotted  is  moved  down  to  blQ  and  the  dotted 

area  is  discarded  area  is  discarded 


Figure  3:  The  three  ways  to  divide  a  cube 
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In  case  (b): 

J3£  =  {(B£\  b')Ub'a}, 

where  bea  =  be; 
W'  =  {(WZ\  w£)Uw£} 

where  wea  =  (w1'*, . . . ,  tu*-1-',  ^'Y',  I/;'*1'', . . . ,  wm''); 
Cf  =  {(<5£\  c')Uc£}. 

In  case  (c): 

where  &£  =  (ft1^, . . . ,  tf-M  ^l±«^i,  &J+M  . . .  ,  &».«); 
W*  =  {(W*\  ^)u<} 

where  u>^  =  u/; 
Ci^{(Ci\  OUcJ,}. 

These  five  steps  constitute  one  complete  iteration  of  the  algorithm.  The  next 
few  lemmata  substantiate  the  claim  that  the  comprehensive  hulls  of  Bn  and  Xn  do 
indeed  give  upper  and  lower  estimates  of  the  production  set.  Lemma  2  says  that 
no  part  of  the  PPF  of  a  firm  that  is  in  some  cube  after  step  3  is  removed  in  step  5. 

Lemma  2.  For  all  ye  6  dYe,    ify  6  UcPC'  C(b,w),   then  y  E  (Jcec  C(b,w). 

Proof/ 

It  must  be  shown  that  no  part  of  the  PPF  of  a  firm  that  is  in  some  cube  after 
step  3  is  removed  as  a  consequence  of  step  5.  Since  only  one  cube  is  altered  by 
step  5,  attention  may  be  focused  there.  Step  5  can  do  three  different  things  to  the 
cube,  c£,  which  is  to  be  divided,  depending  on  the  circumstances.  In  case(a),  cl  is 
removed,  and  two  cubes: 

ci  =  [(6'.', . . . ,  &"•■<);  (w>->, . . .  w'-'->,  h"t+2W"\W'^', ....  »-<)] 

and 

c^[(6'.',,...,6^>.',^±^,6^lf,...,6""');  (»'•',...,»-')] 
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are  added  to  C„  to  form  C„-  But  it  is  clear  that  C(be,we)  =  C(b^,w^)  UC(bL,wi). 
Then  trivially,  since  no  y  £  C(be,we)  is  removed,  no  y  £  dYe  (~)C(b ,  we)  is  removed. 
In  case(b),  it  need  only  be  shown  that  the  parts  of  dYe  inside  the  cube  that  is 
to  be  divided  remain  inside  the  resulting  cube.  Again,  this  is  immediate  since  this 
could  only  be  false  if  there  were  some  ye  6  dYeC\C(be,  w£),  and  yJ,£  <  b  '  ~^w  '  .  But 
this  is  a  violation  of  the  conditions  under  which  (b)  is  invoked.  Similarly,  in  case 
(c),  failure  of  the  Lemma  implies  for  some  ye  6  dYe  (1  C(be,we),  yhl  >  b''  ~^wJ'  , 
which  violates  the  conditions  of  the  case. 


Now  consider  the  following  definitions: 

Jfc  k 

Cn  =  Y,Ci      and      Cn  =  Y,Cen. 
e=i  e=i 

These  objects  are  the  subject  of  the  next  two  corollaries.  Global  cubes  are  not 
actually  collected,  and  do  they  play  any  role  in  the  definition  of  the  procedure,  we 
use  them  briefly  here  to  help  prove  corollaries  2.1  through  2.3.  Corollary  2.1  is  just 
the  global  analog  of  Lemma  2. 

Corollary  2.1  For  all  y  G  dY(c0),  if  y  €  \Jc€Cn  C(b,w)i   thenV  €  \JceCnC(biw)- 

Proof/ 

By  definition,  if: 

ye  |ar(c0)f||  U  C(M) 


then, 


V*  =  l,...,fc,3y<€  \dYef]\    \J  C(b,w)\      s.t.  y  =  ][>' 

c<EC<  J    J  *=1 


21 


Then  by  Lemma  2. 

V^  =  l,...,fc,  yee    |J   C(6,u,), 

and  so  by  definition  of  C„, 

ye    {J  C(6,u;). 

cGCn 


Corollary  2.2  says  that  no  feasible  point  on  the  PPF  is  ever  discarded  if  it  is 
at  least  as  good  as  the  tentative  plan. 

Corollary  2.2  For  all  y  G  dY(c0),  ify  +lo  X  Pn  +  u;,  then  y  G  \JceCn  C{b,w). 

Proof/ 

This  will  be  shown  by  induction.  Since  (JceC  C(b,w)  =  C(bo,w0)  and  po  = 
u>o  =  Yl  wQi  ^or  iteration  0  the  statement  reads:  {y  G  dY(co)  \y  +  u>>:wo+u>}CL 
C(bo,wo).  But  this  is  obviously  true  since 

{y  G  dY(cQ)  |y  +  u)>:u)0+u;}c{y6  C(b0,wQ)  \  y  +  u>  t_  wQ  +  u}  =  C(bo,w0). 

Assume  the  statement  is  true  for  n.  To  show  that  it  is  true  for  n  +  1  we  must 
show  that  the  containment  is  preserved  as  Cn  is  changed  into  Cn+i,  and  as  Cn+i  is 
changed  into  Cn+\-  To  see  this  for  the  first  transition,  take  any  y  G  Ucec  C{b,w) 
and  suppose  that  y  $.  (Jcec  C(b,w).  Then  there  is  a  cube  c  G  Cn  such  that 
y  6  C(b,w)  but  c  ^  Cn+i.  But  then,  according  to  step  2,  if  c  is  removed  from  Cn 
then  6  +  w  -<  pn+i  +  w.  Since  y  <  6  for  any  y  G  dY(co),  it  follows  by  monotonicity 
of  preferences  that  y  +  uj  X  pn+i  +  w.  To  prove  that  containment  is  preserved 
during  the  second  transition,  it  is  sufficient  to  show  that  for  any  y  6  dY(co),  if  y  G 
Ucec  ^^i10),  then  y  G  (Jcec  C(6,tx;).  But  this  is  immediate  from  corollary 
2.1. 
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Corollary  2.3  states  that  the  upper  estimate  at  each  iteration  is  indeed  an 
overestimate  of  all  the  interesting  parts  of  the  PPF. 

Corollary  2.3  For  all  y  G  dY(c0)  ify  +  u>  y  pn  +  w,  then  y  G  UEn. 

Proof/ 

Notice  that  \JceCn  C(b,w)  C  ch(Bn)  =  UEn.  Apply  corollary  2.2. 


Lemma  3  shows  that  the  lower  estimate  at  any  iteration  n  is  as  advertised. 

Lemma  3.  For  all  iterations  n,  LEn  C  Y(co). 

Proof/ 

This  follows  directly  from  the  definition  of  LEn  as  the  comprehensive  hull  of 
feasible  points  and  the  assumption  of  free  disposal. 


Lemma  4  is  a  technical  fact  which  will  be  useful  in  proving  future  lemmata. 
Lemma  4.  For  all  iterations  n,  W„  C  LEln. 

Proof/ 

Assume  that  for  all  n,  and  for  all  ce  G  C£,  there  exists  an  xl  €  X„  such  that 
xl  G  C(be,we).  But  then  for  all  n,  and  for  all  we  G  W^  there  exists  an  xe  G  X„  such 
that  xe  >  we.  Therefore,  Wln  C  ch(Xln)  =  LEen. 

Thus  it  only  remains  to  show  by  induction  that  the  assumption  is  true.  For 
n  =  0,  this  is  obvious  since  Cq  is  the  only  element  of  Cq,  Wq  G  Xq,  and  obviously 
i^o  £  C(6q,Wq).  Now  suppose  that  the  assumption  is  true  for  iteration  n.  Consider 
the  following  two  classes  of  cubes: 
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1)  First  consider  any  cube  ce  £  C„+i  such  that  ce  £  C„.  By  the  induction  hypoth- 
esis, there  is  some  xe  £  Xn  such  that  xi  £  C(be,we).  But  since  be  £  i?£+1  and 
b*  >  xe ,  by  step  3  of  the  algorithm,  xe  £  X*  +  1 

2)  Now  consider  any  cube  ceCn+i  such  that  ce  £"  C„.  This  new  cube  must  have 
been  created  through  the  division  some  cube  ce  £  C„  at  step  5  of  the  algorithm. 
In  case  (a)  of  step  5,  the  firm  reports  a  point  xe  £  /i(c),  and  this  point  is  added 
to  X„+l.  But  the  point  xl  is  in  both  of  the  two  new  cubes  that  result  from 
division  since  it  is  in  on  their  common  boundary.  In  cases  (b)  and  (c),  by  the 
induction  hypothesis,  there  exists  an  xe  £  X„  such  that  xe  £  c.  But  since 
xe  £  dYe,  by  the  hypothesis  of  the  case  x  must  still  be  in  the  one  cube,  c*  that 
results 


4.  Technical  Results  on  Cube  Size 

The  purpose  of  this  section  is  to  give  some  technical  results.  Lemmata  5-8, 
show  that  it  is  possible  to  find  an  upper  bound  on  £^*,  the  length  of  longest  edge 
of  any  cube  associated  with  the  estimate  of  the  Uh  firms  production  set  at  iteration 
n,  and  to  show  that  this  bound  decreases  as  n  goes  to  infinity  in  a  predictable  way. 
The  proofs  of  the  results  are  given  in  the  appendix. 

Recall  that  the  effect  of  fifth  step  of  the  algorithm  is  that  some  cube  is  either 
divided  into  two  cubes  (case  (a)),  or  made  into  a  new  cube  half  as  big  as  the 
original  (cases  (b)  and  (c)).  We  will  call  the  cube  or  cubes  that  result  from  such 
a  division  resultant  cubes.  More  generally,  we  will  want  to  keep  track  of  resultant 
cubes,  cubes  resulting  from  divisions  of  resultant  cubes,  and  so  on.   The  following 
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notation  will  be  used  to  indicate  the  pedigree  of  these  classes  of  cubes.  Consider 
the  set  of  cubes  C„  and  let  some  unspecified  number  of  iterations  pass.  Then  the 
sets  Ci  r,Ce  7Ce  3  . . .  will  refer  to  the  sets  of  cubes  that  are  the  result  of  a  single 
division  of  a  cube  in  C„,  two  successive  divisions,  three  successive  divisions,  etc. 
We  may  now  state  Lemma  5. 

Lemma  5.  For  any  I  =  1, . . . ,  k,  suppose  a  cube  c£  results  from  the  division  of 
a  cube  ct  through  step  5  at  iteration  n  of  the  algorithm.  If  j  =  ILE(ce),  then 
(6*  -  ,^)  =  l(fti  -  ,**)  <  lJ3*\ 

Lemma  6.  For  any  i  =  1, . . . ,  k,  for  any  cfm  €  C„  rm ,  for  all  j  =  1, . . . ,  m,  bJr'm  — 

The  point  of  Lemma  6  is  that  if  a  cube  is  a  result  of  m  divisions  of  some  original 
cube  c* ,  one  of  two  things  must  be  true:  either  each  of  the  m  edges  of  the  resultant 
cube  have  been  divided  exactly  once,  or  at  least  one  of  the  edges  have  been  divided 
twice  (or  more).  In  the  former  case,  all  the  edges  of  the  resultant  cube  are  exactly 
half  the  length  of  those  of  cr.  Thus,  all  the  edges  are  half  as  long  or  less  than  E„*, 
the  longest  edge  of  any  cube  in  C„.  In  the  latter  case,  if  some  edge  is  divided  twice, 
then  the  longest  edge  must  also  have  been  divided  twice  since  it  would  be  the  first 
to  be  divided  a  second  time.  But  then  all  the  edges  of  the  resultant  cube  are  less 
than  or  equal  to  half  the  length  of  longest  edge  of  the  original  cube.  But  this  is  less 
than  or  equal  to  \E^ ,  which  proves  the  Lemma. 

Corollary  6.1  For  any  I  =  1, . . . ,  k,  for  any  ce  6  Ce  m,  for  m'  >  m,  and  for  all 
i  =  l,...,m,  b3'e,  -wJ'e,  <  \El*. 

Corollary  6.1  generalizes  Lemma  6  to  show  that  the  conclusion  holds  for  cubes 
that  are  the  result  of  more  than  m  divisions  of  an  original  cube. 

Lemma  7.  For  any  I  =  1, . . . ,  k,  if  at  some  iteration,  n,  there  are  at  most  Q  cubes 
in  the  set  C„,  then  after  Q(2m  —  1)  more  iterations,  there  will  be  at  most  Q2m 
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cubes  in  the  set  Cr„+o(2»n-i)i  an(^ 

£/[n+Q(2"-l)]   ^   2      n  " 

Lemma  7  extends  the  argument  of  corollary  6.1  to  say  that  if  at  iteration  n, 
there  are  Q  cubes,  then  after  Q(2m  —  1)  more  iterations,  each  of  the  original  Q  cubes 
must  have  been  divided  exactly  m  times,  or  at  least  one  of  the  original  cubes  must 
have  been  divided  more  than  m  times.  In  both  cases  corollary  6.1  may  be  applied 
to  conclude  that  the  longest  edge  of  any  cube  in  Cn+Q(2m-i)  is  at  most  half  the 
length  of  the  longest  edge  of  any  cube  in  C„. 

Lemma  8.  For  any  I  =  1, . . . ,  k,  at  iteration  I{t)  =  (2m  -  1)  £*=0  2am,   EjJt)  < 

i     j?t* 

Lemma  8  builds  on  Lemma  7  to  consider  "blocks"  of  Q(2m  —  1)  iterations  in 
order  to  calculate  how  many  iterations  must  pass  before  E„*  is  smaller  than  ^Eq* . 


5.  Results 

At  last  we  come  to  the  results.  First  we  show  that  the  procedure  is  monotonic 
and  gives  feasible  tentative  plans. 

Theorem  1.   The  set  of  tentative  production  plans,  {p^^-Y  are  feasible,  and  mono- 
tonically  increasing  in  the  preference  order. 

Proof/ 

pn  £  LEn-\  =  ch(Xn-i).  By  monotonicity  of  preferences,  pn  =  x  for  some 
x  £  Xn-i.  But  by  construction,  x  =  ^  xl  for  some  xe  G  Ye  D  Cq  for  £  =  1, . . .  ,  k. 
Thus  pn  =  x  and  so  pen  =  xe  G  Y£Dcq  or  alH  =  1,. . . ,  k.  Thus  for  all  n,  pn  G  Y(c0). 
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Step  1  of  the  algorithm  stipulates  that  pn  £  {x  £  LEn-\  |z+u;>:y  +  <x>VyE 
L^n-i}.  But  Pn-i  is  also  an  element  of  LEn-\  since  only  elements  of  LEn-2  that 
are  strictly  inferior  to  pn-\  are  removed  to  form  LEn-\.  It  is  immediate  from  the 
above  that  pn  >z  pn-\  for  all  n.  Thus,  {pn}^L1  is  monotonically  increasing  in  the 
preference  order. 


There  are  several  ways  to  think  about  convergence.  Theorem  2  shows  that  the 
minimal  distance  between  the  indifference  surface  through  the  tentative  plan  and 
the  optimal  plans  converges  to  zero.  In  addition,  we  show  that  jtqpr~  *s  an  upper 
bound  on  this  distance  at  iteration  n  =  I(t).  This  may  be  taken  as  a  measure  of  an 
upper  bound  on  speed  of  convergence  in  the  preference  ordering  that  is  independent 
of  the  utility  representation.  We  can  use  this  bound  as  a  basis  of  comparison  with 
other  procedures. 

Theorem  2.  Ifn  >  I(t),  u(L(pn  +  w),  U{X*(c0)  +u)  <  ^£ 

Proof/ 

Consider  any  x*  £  X*(cq).  By  corollary  2.2,  x*  £  C(b,w)  for  some  c  £  Cn. 
Then  there  exist  cubes  ce  =  (be,we)  £  Cln  for  £  =  1,...,&  such  that  £^=1  b£  — 
b  >  x*  and  Yle=\  w  ~  w  —  x* ■  But  then  there  are  x*'e  for  £  =  1, . . . ,  k  such  that 
be  >  xim  >  wl  and  ^2(=1  x£*  =  x* .  Since  by  Lemma  4,  we  £  LE„,  it  follows  that: 


!/(££#, *'•)<   ||*"-u>£||   <   ||  6'- 


w 


But  by  Lemma  8,  b*'1  -  w1^  <   ^  for  all  ce  £  C£,  i  =  1, . . . ,  m,  and  £ 
!,...,&.  Then  using  the  definition  of  Euclidean  distance: 


*(£££,**•')< 


\ 


£(*•'  -  "''')2  <  ^/m  (|£) 


Ef* 


ra 

2*+r"' 
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and  so, 

k 


u(LEn,x  )  <  2^     9<+1      -  -^7T1-' 


9<+i 


Butx*  G  J7(X*(x)+u;).  SobyLemmal,i/(L£n,J7(X*(co)+w))  <  ^^.  Similarly, 
by  step  1  of  the  algorithm,  and  Lemma  3,  LEn  +  w  C  L(pn  +  u>).  Thus,  by  Lemma 
1,  u(L(pn  +  «),  U(X*(co)  +  «))  <  §#. 


In  a  planning  environment,  the  CPB  is  usually  assumed  to  know  the  SWF. 
It  is  more  traditional,  as  a  consequence,  to  think  of  convergence  in  utility  terms. 
Theorem  3  shows  that  for  any  utility  representation  of  social  preferences,  the  utility 
of  the  tentative  plans  converges  to  the  utility  of  the  optimal  plans. 

Theorem  3.   Given  any  continuous  utility  representation,  u  :  X  — ►  3?,  of  the  social 
preference  relation,  y,  for  all  x*  G  X*(co),  limn_oo  u(pn  +  to)  =  u(x*  +  u). 

Proof/ 

Assumptions  Al,  and  A2  are  sufficient,  according  to  Debreu(1954),  to  assure 
that  >z  can  be  represented  by  a  continuous  utility  function.  Recall  that  by  Theorem 
2,  for  n  >  I(t),  u(L(pn  +  w),  U(X*(c0)  +  w))  <  Yle=\  2«  +  i '  Thus,  as  n  —>  oo, 
u(L(pn  +  u),U(x*  -f  u)))  —*  0.  Then  we  can  pick  two  sequences  {!„},  and  {Un} 
such  that  for  all  n,  Ln  G  L(pn  +  u),  Un  G  U(pn  +  w),  and  ||  Ln  —  Un  \\  —*  0.  Since 
u  is  continuous,  linin—oo  u(Ln)  —  u(Un)  =  0.  Then,  since  u(Ln)  <  u(pn  +  u)  < 
u(x*  +  u>)  <  u(Un),  limn_^oo  u(pn  +  u>)  -  u(x*  +  u)  =  0. 


Theorem  4  gives  an  upper  bound  on  the  speed  of  convergence  in  utility  terms. 
Such  an  estimate  is  important  because  it  gives  the  CPB  a  basis  to  compare  different 
planning  procedures  and  decide  which  is  best  to  solve  his  specific  problem. 
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Recall  from  real  analysis  that  all  real  continuous  functions  on  a  compact  metric 
space  are  also  uniformly  continuous.  We  will  be  forced  to  strengthen  this  somewhat 
in  order  to  actually  calculate  a  bound  on  the  speed  of  convergence.  In  particular, 
we  will  assume  that  the  SWF  is  proportionally  uniformly  continuous.  In  economic 
terms  this  means  essentially  that  if  two  indifference  curves  are  close  to  each  other 
somewhere,  then  there  is  a  proportionate  bound  on  how  far  apart  they  can  ever  get 
from  one  another  in  the  rest  of  the  goods  space.  This  is  a  bound  of  the  relative 
steepness  and  shallowness  of  the  utility  hill.  This  is  only  a  slight  strengthening  of  the 
existing  assumptions  on  the  economy,  and  as  such  should  be  viewed  as  being  purely 
technical  in  nature.  Given  that  we  already  have  continuity  and  strong  monotonicity, 
we  are  only  eliminating  the  possibility  that  the  utility  hill  has  infliction  points  by 
this  assumption. 
A3)  For  a  given  utility  representation,  u,  of  >:,  3  A  >  0  such  that  Vx  £  C(60,u>o),  ^ 

II  x  —  V  ||     <    ei  then    |  u(x  +  u)  —  u(y  +  u)  \    <  Xe. 

(Proportional  uniform  continuity  of  utility.) 
Since  we  assume  that  the  CPB  knows  his  SWF,  it  is  not  significant  that  he  needs  to 
know  the  A  parameter  of  proportional  uniform  continuity  to  calculate  the  bound  on 
speed  of  convergence  in  utility  terms.  Note  that  A  =  1  if  the  CPB  has  transferable 
utility. 

Theorem  4.  For  any  utility  representation,  u,of>z  satisfying  A3,  any  x*  €  X*(c"o), 
and  any  n  >  I(t),  u(x*  +  u)  —  u(pn  +  u;)  <  A    °+™  . 

Proof/ 

By  Theorem  2  for  any  n  >  I(t),  and  for  all  x*  G  X*(c0),  u(L(pn  +  u>),  U(x*  + 
w))  5-i  2°-M  •  Thus,  for  every  n  >  /(£),  there  exist  xn,  and  yn  such  that  i"  +  uj  E 
U(x* +oj),  yn  +uj  6  L(pn  +uj)  and  ||  xn  —  yn  \\  <  2°(  +  im  •  But  then  by  assumption 
A3,  and  the  fact  the  xn  +  u)  >z  yn  +  w,  u(xn  +  u)  —  u(yn  -f  u)  <  A  -X™  ■  Therefore, 
since  xn+u>>:x*+u  and  yn  +  to  X  pn  +  u>,  the  theorem  is  proved. 
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Finally,  Theorem  5  show  that  the  procedure  converges  not  only  in  utility,  but 
in  addition,  the  tentative  plans  converge  in  quantity  terms  to  actual  optimal  plans. 

Theorem  5.   the  sequence  of  tentative  plans  {pn}  converges  to  the  set  of  optimal 
plans  X*(co). 

Proof/ 

Since  {pn}  is  drawn  from  the  compact  set  C(b0,  io0),  we  need  only  show  that  the 
limit  point  of  every  convergent  subsequence  is  an  element  of  X*(c~o).  So  take  any 
convergent  subsequence  {pt}  and  suppose  that  pt  — *  p* .  For  any  x*  £  X*(c~o),  by 
Theorem  3,  u(pt  +  u)  — ►  u(x*  +u).  Then  by  continuity  u(p*  +  ui)  =  u(x*  +  u>).  By 
the  definition  of  pt,  for  all  t  there  is  a  decomposition  p*,  such  that  for  all  I  =  1, . . . ,  k, 
Pt  G  Y*  nC(&o,iuo)-  Since  Y  nC(&o5  ^o)  is  compact  by  assumption  B2,  there  is  a 
decomposition  p*,  such  that  for  all  I  =  1,. . . ,  k,  pe*  E  Ye  D  C(bo,wo).  Therefore, 
p*  G  Y(co),  and  is  also  undominated  in  the  preference  ordering.  We  conclude  that 
P*  €  X*(cq). 


Unfortunately,  it  does  not  seem  to  be  possible  to  extend  theorem  4  and  find  a 
general  bound  on  the  speed  of  convergence  in  quantity  terms.  To  do  so  we  would 
have  to  know  much  more  about  the  interactions  between  the  preferences  and  the 
feasible  set.  For  some  subclasses  of  economies  (convex  economies  for  example),  it 
may  be  possible  to  find  useful  characterizations  of  these  interactions.  But  this  will 
not  be  attempted  in  the  current  paper. 

Finally,  we  turn  to  the  question  of  stopping  rules.  Except  in  very  special 
cases,  finite  convergence  cannot  be  expected.  So  in  practice,  the  CPB  will  have 
abandon  the  search  at  some  point  and  produce  the  current  tentative  plan.  The 
CPB  must  therefore  devise  a  rule  to  stop  the  procedure  when  then  the  tentative 
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plan  is  sufficiently  "close"  to  an  optimal  plan.  If  the  CPB  is  satisfied  with  defining 
"closeness"  in  utility  terms,  stopping  rules  are  very  easy  to  implement.  All  he 
need  do  is  subtract  the  utility  of  the  tentative  plan  from  the  utility  of  the  most 
preferred  element  of  the  Upper  Estimate,  and  stop  the  search  when  this  number 
falls  below  a  pre-specified  threshold.  If  the  CPB  insists  on  making  a  stopping  rule 
in  quantity  terms,  then  things  are  slightly  more  complicated.  The  CPB  must  find 
an  upper  bound  on  the  distance  in  goods  space  between  the  current  tentative  plan 
and  the  set  of  optimal  plans.  Recall  that  the  algorithm  discards  information  as  it 
progresses.  Then  since  we  know  that  all  X*(co)  E  (Jcec  ^(^iw)i  one  waY  t°  nn(^ 
an  upper  bound  is  to  take  the  maximum  distance  between  the  points  in  this  union. 
Unfortunately,  it  will  not  always  be  the  case  that  this  bound  goes  to  zero  as  the 
number  of  iterations  goes  to  infinity.  So  a  CPB  who  uses  a  quantity  stopping  rule 
like  this  can  never  be  sure  that  he  will  ever  actually  stop.  Stopping  is  more  likely, 
however,  the  closer  set  of  optimal  plans  are  to  one  another.  In  particular,  if  the  set 
of  optimal  plans  is  a  singleton,  stopping  is  always  guaranteed. 


6.  Conclusions 

In  conclusion,  the  new  procedure  is  well  defined,  monotonic,  convergent,  and 
even  in  a  finite  number  of  iterations,  if  the  true  production  set  is  a  step  function.11 
In  addition,  The  message  space  is  simple,  and  the  large  class  of  comprehensive 
and  monotonic  economies  can  be  treated.  An  important  feature  of  the  procedure 
presented  in  this  paper  is  its  ability  to  discriminate  against  irrelevant  parts  of  the 


In  Malinvaud's  and  Weitzman's  procedure,  finite  convergence  is  obtained  when  the  production  set 
in  polyhedral.  This  is  because  in  this  case  the  production  set  can  be  exactly  approximated  by  the 
estimates  that  the  procedure  generates.  The  same  thing  holds  here  for  step  function  production 
sets  which  can  be  exactly  approximated  by  Quantity-Quantity  type  procedures.  This  is  not  a  very 
important  case,  however,  and  no  proof  will  be  offered  of  this  assertion. 
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goods  space  and  so  narrow  the  area  of  search.  It  is  unfortunate  that  there  is  no 
obvious  way  to  incorporate  this  into  the  estimate  of  the  speed  of  convergence.  The 
greatest  advantage  of  this  procedure  is  that  its  design  makes  it  makes  it  possible 
to  directly  calculate  an  upper  bound  on  the  number  of  iterations  it  takes  for  the 
tentative  plans  to  be  within  any  given  e  in  utility  of  the  optimal  plans.  Although 
this  is  a  very  loose  upper  bound,  it  still  can  serve  as  a  basis  of  comparison  for  speed 
of  convergence  with  other  procedures.  This  is  essential  if  a  CPB  is  to  make  an 
informed  decision  about  which  planning  procedure  to  use.  In  addition,  at  any  given 
iteration,  it  is  possible  to  place  an  upper  bound  on  how  far  the  current  tentative 
plan  is  from  an  optimal  plan  both  in  utility  terms,  and  in  the  actual  goods  space. 
This  allows  the  CPB  to  make  a  reasonable  judgment  about  whether  the  current 
plan  is  good  enough,  or  to  continue  searching. 


Appendix 

Lemma  5.  For  any  £  =  l,...,k,  suppose  a  cube  c\.  results  from  the  division  of 
a  cube  ce  through  step  5  at  iteration  n  of  the  algorithm.  If  j  =  ILE(ce),  then 
(#t  -  wi/)  =  I(6i  _  WJ)  <  lEi\ 

Proof/ 

Consider  any  cube  cf  that  results  from  a  division  along  the  jth  edge  of  a  cube 
ce.  The  cube  c(r  must  have  resulted  from  application  of  one  of  cases  (a),  (b),  or  (c) 
in  step  5.  In  case  (a),  cer  can  take  two  forms: 

(*'•',...>' V*),  («,'•',...,  ^±^,...,  **■«) 

or 

'(^ '"^ ^.(^ ^.-.'^ 

In  case  (b),  only  one  form  is  possible: 
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(6     ,...  ,6-/'  ,...,6       ),{w      , . 
Likewise,  in  case  (c),  cr  must  take  the  form: 


bi>e  +  w^e 


1£  V<e  +  wM  t         ie 


Then  since  either, 
(i)  bp  -  wpe  =  b^  -  ^|u±l  =  bLi^eLL  or 


WJ' 


bi.(-wJ-{ 


w 


i,e  _  6'-<W-< 


10 


m,^ 


U,m'') 


(ii)  H'e  -  *>*/  = 

But  since  a  cube  is  always  divided  along  it  longest  edge,  and  by  definition,  no  edge  of 
any  cube  can  be  longer  that  E^  at  iteration  n,  (bJr'e  —  wJr'e)  =  UbJ'e  —  w3'e)  <  \El* . 


Lemma  6.  For  all  £  =  1, . . . ,  k,  for  any  c£m  €  Cen  rm ,  for  all  j  =  1, . . .  ,  m,  bJr'n 


w 


j,t 


<  \E<n\ 


Proof/ 


We  start  by  distinguishing  two  exhaustive  subclasses  of  cubes  in  C„  rm. 

1.  First  consider  cubes  cerTn  E  C„  rm  for  which  no  edge  has  ever  been  subjected 
to  two  separate  divisions.  That  is,  cubes  that  are  the  result  of  m  divisions 
of  an  original  cube  in  C£,  and  for  which  each  of  the  m  edges  has  been 
divided  exactly  once.  By  Lemma  5,  (b3r'm  —  wJr'm)  =  ^{b3,e  —  whl)  for  all 
j  =  1, . . . ,  m.  But  by  definition,  (b3,1  —  w3,i)  <  E„*  for  all  j  =  1, . . . ,  m. 
Thus,  for  cubes  in  this  class,  (bJr'm  —  wJr'm )  <  \E^  for  j  =  1, . . . ,  m. 

2.  Now  consider  cubes  cfm  €  Cln  rm  for  which  at  least  one  edge  has  been 
subjected  to  at  least  two  separate  divisions.  Without  loss  of  generality, 
suppose  that  edge  j  is  the  longest  edge  of  the  original  cube  c^,  and  so  is 
the  first  to  be  divided.  Then  by  Lemma  5,  (b3r,e  —  w3r'e)  =  ^{b3'e  —  w3,e). 
But  by  assumption,  ^{b3'*  -  w3<e)  >  \{b{>1  -  w1*1)  for  all  i  jk  j.  So  if 
any  edge  is  divided  twice  through  step  5,  then  the  ]th  edge  is  also  divided 
twice.  But  since  only  the  largest  edge  is  ever  subject  to  division,  \E^  > 
i(6J'^  —  wJ'e)  >  (6j.'m  —  w^m )  for  all  i  =  1, . . . ,  m.  Thus,  for  cubes  in  this 
class  as  well,  (bJr'm  —wJr'm )  <  hE„*  for  j  =  1, . . . ,  m.  The  Lemma  is  proven. 


Corollary  6.1  For  all  I  =  1, . . . ,  k,  for  any  ce  (E  Ce      ,  for  m'  >  m,  and  for  all 


j  =  l,...,m,  o^, 
Proof/ 


w 


i,l 


n,rm 
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By  Lemma  6,  for  any  c£m  €  C„  rm,  for  all  j  =  1, ...  ,m,  6^'m  —  u;£in  <  |-E"n*. 
Since  it  is  impossible  for  any  edge  of  any  cube  to  be  increased  as  a  result  of  a 
division,  b3'  ,  —  w3'  ,  <  \Ei*  for  m'  >  m. 


Lemma  7.  For  all  £  =  1, . . . ,  k,  if  at  some  iteration  n  there  are  at  most  Q  cubes  in 
the  set  C„,  then  after  Q(2m  —  1)  more  iterations,  there  will  be  at  most  Q2m  cubes 
in  the  set  C*n+Q(2m_l^  and 

£/[n+Q(2--l)]   ^    2      n  ' 

Proof/ 

The  first  part  is  easy  to  show.  Each  iteration  can  add  at  most  one  cube.  This 
is  because  the  step  3  eliminates  cubes,  if  it  does  anything  at  all,  and  (a)  of  step  5 
adds  one  additional  cube  while  (b)  and  (c)  leave  the  number  of  cubes  unchanged. 
Thus,  an  upper  bound  on  the  number  of  cubes  at  the  end  of  iteration  n  +  Q(2m  —  1) 
when  there  were  at  most  Q  cubes  at  iteration  n  is  Q  +  Q(2m  —  1)  =  Q2m . 

To  see  the  second  part,  suppose  initially  that  each  of  the  cubes  in  C„  is  divided 
once  before  any  of  the  cubes  is  redivided.  Then  after  Q  iterations,  C„  D  C?.q,  =  0. 

Also,  Cf+Q,  will  consist  of  at  most  2Q  cubes.    Now  let  each  cube  in  C*+q  be 

divided  once  before  any  is  divided  a  second  time.  Since  C?n  ,  q,  contains  at  most 
2Q  cubes,  this  will  take  at  most  2Q  more  iterations.  At  the  end  of  these  iterations 
Cf  +Q+2Q]  w*^  contain  at  most  Q22  cubes,  and  by  construction,  CnnC[n+Q+2Q]  =  0 
and  C„  r  Pi  Cf  ,  Q+2Q1  =  ®-  Suppose  that  this  process  continues,  and  each  cube 
in  the  set  C„  rX  is  divided  before  any  cube  in  C^  rX  +  l  is.  Then  by  the  end  of 
n  +  Q  +  2Q  +  Q22  +  Q23  +  . . .  +  Q2m~l  iterations,  there  are  at  most  Q2m  cubes 
and 

for  all  x  <  m.  Note  that 

0(1  +  2  +  22+, . . . ,  +2™-1 )  =  Q(2m  -  1). 
We  conclude  that  all  elements  of  Cf  ,n<nm    ^^^  are  elements  of  some  Ce      .  where 

ln  +  y(^™  — 1)1  n,rm 

m'  >  m.  But  by  corollary  6.1,  for  any  c£  €  Ce  m,  for  all  m'  >  m,  and  for  all 
j  =  1, . . . ,  m,  b3'e,  -  w3'e,  <\El*.  Therefore: 


£'[n  +  Q(2"»-l)]   ^   2-^"   " 
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Now  suppose  that  things  do  not  develop  so  neatly  and  some  cubes  in  C^  rX+1 
are  divided  before  some  in  Cln  rX .  Then  assume  some  cubes  have  been  divided  less 
than  m  times  by  the  end  of  iteration  n  +  Q(2m  —  1)  (or  else  corollary  6.1  can  be 
applied  directly  as  in  the  first  case).  Then  since  one  cube  must  be  divided  at  each 
iteration,  if  some  cube  is  divided  less  than  m  times,  some  other  cube  must  have 
been  divided  more  than  m  times  by  the  end  of  iteration  n  +  Q(2m  —  1).  Thus  at 
some  stage,  (say  iteration  n'  where  n  <  n'  <  n  +  Q(2m  —  1)),  some  ce  (E  Cln  rm 
is  divided.  But  this  can  only  happen  if  for  all  cl  6  C*,,  and  all  i  =  l,...,m, 
bl,i  —  wl'£  <  bi'e  —  w3,e,  where  j  =  ILE(ce).  Since  further  divisions  of  cubes  after 
iteration  n'  will  not  increase  the  length  of  any  edge  of  any  cube,  and  by  Corollary 
6.1  &>•'  -  wJ'e  <  ±Een*,  we  conclude: 


^[n+Q(2--l)]    S   ^^n  • 


Lemma  8.  For  all  I  =  1,. .  .,k,  at  iteration  I(t)  =  (2m  -  l)£l=02am,   Ejm(t)  < 


i    jpe* 


2t  +  i  "^0    • 

Proof/ 

This  is  shown  by  induction.  For  t  =  0,  I(t)  =  2m  —  1,  so  it  must  be  shown  that 
Efcm_l]  <  2^0*-  At  iteration  n  =  0,  there  are  Q  =  1  cubes  in  Cg,  so  by  Lemma  7, 
after  1  x  (2m  —  1)  additional  iterations, 

j?e*  —  j?**  <^  i  j??* 

■C/[0  +  (2m-l)]   -  rj[(2m-l)]   -    l111^    • 

Now  assume  that  the  statement  is  true  for  t.  Then  EjJt^  <  ^tTT-^o*-  Since  each 
iteration  adds  at  most  one  cube,  and  there  exists  only  one  cube  at  iteration  0,  after 
I(t)  iterations,  there  are  at  most 

l{t)  +  1  =  (2m  -  1)  (  J2  23m  )  +  l  =  (  £  rm  ~  J2  29m  )  +  l  =  2m(<+1) 

\s=0  /  \s=l  3=0  / 

cubes  in  Cj,ty  Then  by  Lemma  7,  after  at  most  2m(*+1)(2m  —  1)  more  iterations, 

r/[/(0+2mC  +  1)(2"'-l)]  -  2^/(0  -  2        2t+1     °  2*  +  2     ° 

But 

I(t)  +  2m(<+1)(2m  -  1)  =  (2m  -  1)  [  J2  rm  )  +  (2m  -  l)2m(<+1)  = 
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(2m-l)(^2sm  J  =  I(t  +  l). 
So  at  I(t+1),  Efit+1)  <  ifeEZ*. 
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